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Abstract: Many desert expressways are affected by the deposition of the wind-blown sand, which might block 
the movement of vehicles or cause accidents. W-beam central guardrails, which are used to improve the safety 
of desert expressways, are thought to influence the deposition of the wind-blown sand, but this has yet not to 
be studied adequately. To address this issue, we conducted a wind tunnel test to simulate and explore how the 
W-beam central guardrails affect the airflow, the wind-blown sand flux and the deposition of the wind-blown 
sand on desert expressways in sandy regions. The subgrade model is 3.5 cm high and 80.0 cm wide, with a 
bank slope ratio of 1:3. The W-beam central guardrails model is 3.7 cm high, which included a 1.4-cm-high 
W-beam and a 2.3-cm-high stand column. The wind velocity was measured by using pitot-static tubes placed 
at nine different heights (1, 2, 3, 5, 7, 10, 15, 30 and 50 cm) above the floor of the chamber. The vertical 
distribution of the wind-blown sand flux in the wind tunnel was measured by using the sand sampler, which 
was sectioned into 20 intervals. In addition, we measured the wind-blown sand flux in the field at K50 of the 
Bachu-Shache desert expressway in the Taklimakan Desert on 11 May 2016, by using a customized 
78-cm-high gradient sand sampler for the sand flux structure test. Obstruction by the subgrade leads to the 
formation of two weak wind zones located at the foot of the windward slope and at the leeward slope of the 
subgrade, and the wind velocity on the leeward side weakens significantly. The W-beam central guardrails 
decrease the leeward wind velocity, whereas the velocity increases through the bottom gaps and over the top 
of the W-beam central guardrails. The vertical distribution of the wind-blown sand flux measured by wind 
tunnel follows neither a power-law nor an exponential function when affected by either the subgrade or the 
W-beam central guardrails. At 0.0H and 0.5H (where H=3.5 cm, which is the height of the subgrade), the sand 
transport is less at the 3 cm height from the subgrade surface than at the 1 and 5 cm heights as a result of 
obstruction by the W-beam central guardrails, and the maximum sand transportation occurs at the 5 cm height 
affected by the subgrade surface. The average saltation height in the presence of the W-beam central guardrails 
is greater than the subgrade height. The field test shows that the sand deposits on the overtaking lane leeward 
of the W-beam central guardrails and that the thickness of the deposited sand is determined by the difference 
in the sand mass transported between the inlet and outlet points, which is consistent with the position of the 
minimum wind velocity in the wind tunnel test. The results of this study could help us to understand the 
hazards of the wind-blown sand onto subgrade with the W-beam central guardrails. 
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1 Introduction 


The W-beam central guardrails are commonly used on desert expressways to contribute to vehicle 
safety (Jason et al., 2005). In recent years, the deposition of the wind-blown sand on the subgrade 
surface with the W-beam central guardrails has become a concern (Li et al., 2016) because 
sand-blown erosion and sand deposition are the primary sources of sand damage on desert 
expressways (Lei et al., 2008; Cheng et al., 2015). The height and slope of the subgrade affect the 
airflow, which correlates negatively with wind velocity, and eddy airflow forms leeward of the 
subgrade (Yang et al., 2010). The settlement of the sand particles at the foot of the leeward and 
windward slopes results in sand deposition on the roadbed (Zhang et al., 2010; Yang et al., 2011). 

Taklimakan Desert, the largest mobile desert in China, is characterized by arid climate with little 
precipitation. The Taklimakan Desert Highway which goes across the Taklimakan Desert was built in 
1995. Along this highway, a comprehensive sand-control system of straw checkerboard barriers and 
fences was installed, which can decrease the wind velocity and sand deposition, so as to protect the 
highway from being buried by the wind-blown sand (Cheng et al., 2015). Huang et al. (2013) and Xu 
et al. (2018) also simulated the turbulent flow and the motion of the sand particles over the surface of 
straw checkerboard barriers. Wind fences, which are used to protect areas from high wind action and 
sand damage, have been studied to determine the distance of optimum protection, the wind reduction 
efficiency, and the optimum porosity (Dong et al., 2006, 2007; Li and Douglas, 2015; Zhan et al., 
2017). The shelterbelt systems in the Taklimakan Desert have also been studied, and the results show 
that the shelterbelt can not only control the wind-blown sand damage but also provide ecological 
benefits (Lei et al., 2008). 

Significant research has been done to understand how fences affect the vertical distribution of the 
wind-blown sand flux, how fence porosity affects the sediment deposition and why sand 
accumulates upwind of low-porosity fences but downwind of high-porosity fences (Dong et al., 
2004; Cornelis and Gabriels, 2005). Large-diameter sand particles deposit on the windward side of 
fences, whereas suspended sand particles pass through the fences, and the sand particles deposited on 
the leeward side of windbreaks stay put because of the direction and intensity of the eddy airflow 
(Zheng et al., 2011). The mean sand velocity decreases dramatically and the particle concentration 
decreases leeward of a fence, further, the vertical distribution of the sand around fences follows an 
exponential function (Zhang et al., 2010). 

In desert regions, with the development of road transportation, there is a growing demand for the 
building of desert expressways (Zheng et al., 2011). High-speed expressways differ substantially 
from low-speed roads in terms of subgrade structure and safety requirements. Higher subgrade, 
central guardrails, and the anti-glare nets of desert expressways tend to increase the sand deposition 
on the subgrade surface (Yang et al., 2010). The W-beam central guardrails have generated new types 
of sand damage, but their interaction with airflow and wind-blown sand remains unclear. This study 
uses wind tunnel test and field test to investigate the airflow associated with the wind-blown sand 
under the W-beam central guardrails. The results of this study could help us to understand the 
wind-blown sand hazards caused by the subgrade with the W-beam central guardrails. 


2 Methods 


2.1 Wind tunnel test of wind velocity 


A scaled simulation experiment was conducted in a wind tunnel at the Desert Research Station of 
Mosuowan, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences. The total 
length of the blower-type non-circulating wind tunnel was 16.2 m, which included an 8-m-long test 
section. The test section was 1.3 m wide and 1.0 m high. The inlet wind velocity at the wind tunnel 
ranged from 1 to 20 m/s and the boundary layer in the case of pure wind was about 15 cm for an inlet 
velocity of 12 m/s (Zheng et al., 2012). 

Figure 1 shows the experimental setup. The W-beam central guardrails combined a 28-cm-high 
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W-beam and a 46-cm-high stand column, as used on the Bachu-Shache desert expressway in the 
Taklimakan Desert. The wind tunnel model was a 1:20 scaled-down version of a prototype 
expressway design. The model consisted of two parts: the subgrade and the W-beam central 
guardrails. The subgrade model was 3.5 cm high and 80.0 cm wide, with a bank slope ratio of 1:3. 
The W-beam central guardrails model was 3.7 cm high, which included a 1.4-cm-high W-beam and a 
2.3-cm-high stand column. Figure 2a shows a photograph of the W-beam central guardrails in the 
field and Figure 2b shows the guardrails used in the wind tunnel test. The wind velocity was 
measured in two cases: (1) subgrade without the W-beam central guardrails and (2) subgrade with the 
W-beam central guardrails. Specifically, the wind velocity was measured by using pitot-static tubes 
placed at nine different heights (1, 2, 3, 5, 7, 10, 15, 30 and 50 cm) above the floor of the chamber. 
The pitot-static tubes were connected to calibrated pressure transducers via thin plastic tubing. The 
test points in front of the W-beam central guardrails were —20.0H, —15.0H, —10.0H, —8.5H, —7.0H, 
3.5H and —1.0H (where H=3.5 cm, which is the height of the subgrade), the test points in lee of the 
W-beam central guardrails were 1.0H, 3.5H, 7.0H, 8.5H, 10.0H, 15.0H, 20.0H, 25.0H and 30.0H, 
and the test point in central of the W-beam guardrails was 0.0H. The inlet wind velocity was 
determined by hot-wire anemometry at the height of 50 cm and with the wind velocities of 8, 10 and 
12 m/s in the wind tunnel test. The temperature, relative humidity, and atmospheric pressure were 
measured by using a hygrothermograph just outside the wind tunnel. 


Hot-wire anemometry 450 


Pressure sensor 


Data acquisition system 
Fig. 1 Schematic diagram of experimental setup 


The physical wind flow in the wind tunnel test can be understood by applying the correct 
simulation principles. The physical processes of the full-scale flow can be simulated in the 
laboratory by matching the non-dimensional coefficients of the equations of motion in the scaled 
laboratory experiment to those obtained in the field wind tunnel (White, 1996). Three parameters 
need to be considered: the Rossby number, the bulk Richardson number and the Reynolds number 
(Re=pvL/u, where p is the density of the fluid (cm?), v is the velocity (m/s), L is the characteristic 
length (m) and yw is the viscosity coefficient (Pa-s)). The effects of the Rossby number can be 
ignored if the laboratory test is restricted to the lowest 10.0%—15.0% of the height of the boundary 
layer. If the laboratory test is done in the lowest 7.2% of the height of the boundary layer, then the 
longitudinal velocity spectrum in the inertial subrange is modeled accurately. The bulk Richardson 
number is zero for neutral-stability atmospheric flow, which is easily simulated in the laboratory by 
isothermal flow. A Reynolds number greater than 1.0x10* is independent of scale. In this study, we 
used a minimum wind velocity of 8 m/s in the wind tunnel test, with a subgrade width of 80 cm. 
This gave a Reynolds number of 46.0x10*. The Reynolds number was therefore independent so the 
airflow over the scale model was dynamically similar to that over the full-scale field versions 
(White, 1996). We also considered the geometric similarity between the prototype and the scale 
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model and matched the shape and materials of the model to those of the prototype. We set the 
degree of congestion of the model size being less than 6.0% to avoid interference by the wind 
tunnel wall. The maximum windward area of the tested model was 392.0 cm? and the degree of 
blocking was 3.9%. 


(a) 


Fig. 2 Photograph of the W-beam central guardrails on the Bachu-Shache desert expressway in the Taklimakan 
Desert (a) and overview of the scaled-down W-beam central guardrails in the wind tunnel test (b) 


2.2 Structure of the wind-blown sand flux and patterns of the sand deposition 


The structure of the wind-blown sand flux was measured in the wind tunnel. The median diameter of 
the sand particles was 138.8 um, as measured by a Malvern particle-size analyzer (Mastersizer 2000, 
Malvern Instruments Ltd., Malvern, UK). The sand particles were laid (in a 20 cm deep layer and 
about 1 m long and 1 m wide) on the wind tunnel floor about 2.0 m downwind from the entrance of 
the test section in the wind tunnel. The sand sampler was used to measure the vertical distribution of 
the sand particles and the wind-blown sand flux. The sampler was 20 cm tall and sectioned into 20 
intervals (area of each interval: 1 cmx1 cm), and its shape and working principles were the same as 
reported by Dong et al. (2003). The sand sampler was positioned at —20.0H, —10.0H, —8.5H, —7.0H, 
0.0H, 7.0H, 8.5H, 10.0H and 20.0H in the wind tunnel to measure the properties of the subgrade and 
at —7.0H, —0.5H, 0.0H, 0.5H and 7.0H to determine the effect of the W-beam central guardrails. The 
mass of the sand particles collected in the sand sampler was measured by an electronic balance. The 
experimental wind velocity was 8 m/s and the test lasted for 10 min. 

The field test for wind-blown sand flux was conducted at K50 of the Bachu-Shache desert 
expressway on 11 May 2016. This desert expressway was 28 m wide and had the W-beam central 
guardrails. The wind-blown sand flux was measured by a customized gradient sand sampler, which 
was designed by Wang et al. (2016). The sand sampler had 14 gradients (at heights of 2.0, 5.0, 7.5, 
9.0, 13.0, 20.0, 25.0, 30.0, 37.0, 40.0, 50.0, 60.0, 70.0 and 78.0 cm) and a sand collector with the 
cross-sectional area of 1.7 cmx1.7 cm. The sand sampler was used to make 10-min measurements on 
the sand bed, the windward shoulder, and the leeward shoulder of the expressway to test the 
transportation of the wind-blown sand in the field. The thickness and area of the sand deposits after 
the sandstorm were measured by using a flexible ruler. 


2.3 Data analysis 


2.3.1 Reduction rate of the wind velocity 


The wind velocity field is a comprehensive index indicating the movement of the airflow. The wind 
velocity at a given height and distance above the subgrade may be expressed in terms of the reduction 
rate in the horizontal wind velocity Rc(x,z) (Wang et al., 2017): 

Re... 2) = A= tyz) Í Ug x,2)) * 100% P (1) 


where x is the horizontal distance from the guardrails (H); z is the height above the surface (cm); u¢,) 
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is the wind velocity above a surface with either subgrade or the wind velocity at the measuring point 
(m/s); and woz) is the inlet wind velocity (m/s). Here Rcq)<0 indicates that the wind velocity is 
higher than the inlet velocity, and Rc&,>0 indicates that the wind velocity is lower than the inlet 
velocity. 

2.3.2 Energy loss calculation 


The energy loss of the airflow (AZ) in a unit volume can be expressed by Equation 2 (Ni and Li, 
2001): 


AE= PCs T Hea k (2) 


where AZ is the energy loss of the airflow (J/m?); p is the air density (kg/m°); wow.) is the inlet wind 
velocity (m/s); and uz) is the wind velocity at the measuring point (m/s). A larger AE corresponds to 
less wind power, and AF<0 indicates increasing wind power. 

2.3.3 Calculation of the sand deposit thickness 

For the desert expressway, the windward side of the road shoulder acts as an inlet and the leeward 
side of the shoulder acts as an outlet. The difference in mass of the wind-blown sand transport 
between the windward and leeward sides of the expressway gives the mass of the sand deposited on 
the subgrade surface. If the area of the sand deposition on the subgrade is 1.0 m wide, then the 
thickness of the sand deposition layer can be determined by: 


d=- 2) Z0) xt, (3) 


p x10000 
where d is the thickness of the sand deposition (cm); Q1 is the mass of the sand entering the system 
from the windward shoulder (g/(h-m)); Q2 is the mass of the sand leaving the system from the leeward 
shoulder (g/(h-m)); p is the bulk density of the sand particles (g/cm+); and ż is the duration during 
which the wind-blown sand accumulated (h). 
2.3.4 Statistic analysis 


We used Origin Pro 9.0 for the linear regression analysis and to plot the results. Surfer 12.0 was used 
to draw the wind velocity field and AutoCAD version 2007 was used to draw the schematic diagram 
of the experimental setup. 


3 Results 


3.1 Wind velocity field 


Figure 3 shows the horizontal iso-velocity field of the subgrade with and without the W-beam 
guardrails, which is measured under a free stream at wind velocities of 8, 10 and 12 m/s. For the 
subgrade without the W-beam central guardrails, the changes in the horizontal wind velocity can be 
divided into five regions (Fig. 3a). The first region (labeled I in Fig. 3a) is a region of decreasing 
horizontal wind velocity in front of the subgrade. A slight reduction occurs in the horizontal wind 
velocity over the subgrade in the range from —20.0H to —15.0H, where the reduction coefficient of the 
horizontal wind velocity is <9.0%. As the distance from the subgrade surface decreases, the 
horizontal wind velocity decreases at the height of 1 cm above the subgrade surface, with the 
reduction coefficients between 28.0% and 31.0% in the distance range from —12.0H to —11.0H. 

The horizontal wind velocity increases in the second region (labeled II in Fig. 3a), where the 
airflow is compressed along the windward slope and the wind velocity increases by about 3.0%. The 
maximum wind velocity occurs on the windward shoulder of the subgrade. The airflow is diffused 
and passes smoothly over the subgrade surface, and the horizontal wind velocity is stable. The third 
region (labeled III in Fig. 3a) is a region of stable horizontal wind velocity over the subgrade surface. 
The reduction coefficients of the horizontal wind velocity are 0.2%-—0.6% at heights greater than 15 
cm above the subgrade surface and 10.6%-—18.0% at heights less than 15 cm above the subgrade 
surface. The decrease in the horizontal wind velocity near the subgrade surface is a result of the 
surface drag force. The fourth region (labeled IV in Fig. 3a) is the one of decreasing horizontal wind 
velocity in the lee of the subgrade, where airflow separation occurs and a low-velocity region forms 
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from 10.0H to 11.5H. The horizontal wind velocity decreases to a minimum of zero at the foot of the 
leeward slope and the reduction coefficient decreases with the height above the subgrade surface. 


Wind velocity (m/s) 


Height (H) 


10.0 


Height (H) 


Height (A) 


20.0 15.0 10.0 5.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 
8.0 
6.0 
4.0 
2.0 
0.0 


Height (H) 


—20.0 -15.0 i X; 5.0 10.0 i 20.0 25.0 30.0 
Distance (H) 


Fig. 3 Horizontal iso-velocity contours of the subgrade without the W-beam central guardrails at the inlet wind 
velocity of 12 m/s (a), and horizontal iso-velocity contours of the subgrade with the W-beam central guardrails at the 
inlet wind velocities of 12 m/s (b), 10 m/s (c) and 8 m/s (d). H is the height of the subgrade (H=3.5 cm). I, I, HI, ITV 
and V represent the five regions of the horizontal wind velocity affected by the subgrade. 


In the fifth region (labeled V in Fig. 3a), the wind velocity recovers. Specifically, the wind velocity 
increases downwind of the point of reattachment of the separated flows, and the horizontal wind 
velocity recovers with increasing distance from the subgrade surface. The subgrade affects the 
airflow up to 20.0H. After this point, the reduction coefficient of the horizontal wind velocity is less 
than 10.0% and the airflow gradually recovers to its inlet velocity. 

The horizontal iso-velocity contours of the subgrade with the W-beam central guardrails at the 
inlet wind velocities of 12, 10 and 8 m/s are shown in Figures 3b—d. Under the obstruction of the 
subgrade with the W-beam central guardrails, the wind velocity decreases in the range from —3.5H to 
30.0H, and the reduction coefficient of the horizontal wind velocity reaches 47.0% at the height of 
3.5 cm above the subgrade in 0.0H. The reduction coefficients are 27.0% at the test point —3.5H 
above the subgrade surface and 18.0% at the test point —0.5H at the height of 1.5 cm above the 
subgrade surface. The reduction coefficient increases to 47.0% at the height of 3.5 cm in front of the 
W-beam central guardrails and the horizontal wind velocity increases by 1.0% to the inlet velocity at 
the height of 6.5 cm above the W-beam central guardrails. The reduction coefficient is 37.0% at the 
test point 0.5H behind the W-beam central guardrails at the height of 1.5 cm above the subgrade 
surface (Fig. 3b). The minimum horizontal wind velocity on the subgrade surface varies in the range 
of 3.5H—3.7H, decelerating to 0.93, 1.16 and 1.41 m/s at the inlet wind velocities of 8, 10 and 12 m/s, 
respectively; and the reduction coefficient reaches 60.0% at the height of 1.5 cm above the subgrade 
surface (Figs. 3b—d). At the 30.0H test point, the horizontal wind velocity recovers to 72.0%—80.0% 
of the inlet velocity. As the inlet wind velocity increases, the horizontal wind velocity-recovery 
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distance decreases. 
3.2 Changes in the energy loss of the airflow 


Under the influence of the subgrade, the energy loss of the airflow varies from —16.90 to 60.45 J/m? 
(Fig. 4). The increase in the airflow energy is mainly concentrated at the height of 5 cm above the 
expressway shoulder, where wind erosion is easy to happen because of the increase in the airflow 
energy. Under the influence of the subgrade and W-beam central guardrail structure, the energy loss 
of the airflow varies from —2.90 to 90.40 J/m’, and the airflow energy increases above and below the 
W-beam central guardrails. At 0.5H behind the W-beam central guardrails, the energy loss of the 
airflow is 90.40 J/m?, so the wind-blown sand particles are easily unloaded at this point. Overall, the 
energy loss of the airflow is greater for the subgrade with the W-beam central guardrails than for the 
subgrade without the W-beam central guardrails (Fig. 4), and the power of the airflow decreases so 
that the wind-blown sand particles settle in this region. 
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Fig. 4 Energy loss of the airflow for the subgrade with the W-beam central guardrails and without the W-beam 
central guardrails under different distances from the subgrade at the vertical heights of 1 cm (a), 3 cm (b), 5 cm (c), 
7 cm (d), 10 cm (e) and 15 cm (f) 


3.3 Changes in vertical distribution of the wind-blown sand flux 


In the presence of the subgrade and W-beam central guardrails, the vertical distribution of the 
wind-blown sand flux does not decay exponentially with height, and the maximum sand transport is 
concentrated higher above the subgrade (Fig. 5). The subgrade (without the W-beam central 
guardrails) causes the maximum sand transport to concentrate between 3 and 5 cm heights above the 
subgrade surface and to decrease in the heights of 0—1 cm (Fig. 5a). The characteristics of the sand 
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transport can be attributed to the underlying solid surface of the subgrade. For this particular 
subgrade surface, the saltating particles lose only small amounts of momentum in collisions with the 
solid subgrade surface and the majority of the saltating particles jump higher than the level of the soft 
sand bed. At the toe of the leeward slope, the sand particles are lifted by the subgrade and the 
maximum transport of the sand occurs at the 5 cm height above the subgrade surface. At the front and 
on the leeward side of the subgrade, the sand transport decreases with increasing height, which 
conforms to the structure of the sand flux on a flat sand bed (Fig. 5a). 

As shown in Figure 5b, the W-beam central guardrails change the structure of the wind-blown 
sand flux over the expressway subgrade. Under the influence of the W-beam central guardrails, the 
transport of the sand decreases from 65.0% (—7.0H) to 30.0% (0.5H) over a height of 1 cm above the 
subgrade surface. Some sand particles are lifted by the up-flow, and the transport of the sand 
increases with height. However, although some of the sand particles are lifted, the maximum 
transport of the sand is concentrated at the height of 1 cm above the subgrade surface as a result of the 
bleed flow. At test points of 0.0H and 0.5H, the transport of the sand at the height of 3 cm above the 
subgrade surface is less than that at the heights of 1 and 5 cm above the subgrade surface, 


respectively, which is attributed to the obstruction by the W-beam central guardrails. 
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Fig. 5 Structure of the wind-blown sand flux at different test points as a function of height above the subgrade 
surface. (a), subgrade without the W-beam central guardrails; (b), subgrade with the W-beam central guardrails. 


The average saltation height characterizes the profile of the wind-blown sand flux, which is 
defined as the height corresponding to the 50% cumulative percentile (Wu, 1987). Figure 6 shows the 
average saltation height on the subgrade (without the W-beam central guardrails) surface and on the 
subgrade with the W-beam central guardrails. The average saltation height is 0.61—0.89 cm in front 
of the subgrade and the windward shoulder, and then increases to 2.72-3.96 cm on the subgrade 
surface of 0.0—-7.0H. The average saltation height is 4.31-5.17 cm in the leeward of the subgrade. 
The saltation flux over the subgrade surface is maintained by a successive series of saltating 
particles; they obtain their energy from the wind while they are in moving air and then collide 
elastically with the subgrade and rebound with more energy from the air so that they rise to a greater 
height. The average saltation height for the subgrade with the W-beam central guardrails is greater 
than that for the subgrade without the W-beam central guardrails. It varies from 3.94 to 5.06 cm, 

which suggests that more particles move in a higher layer near the guardrails. This result is attributed 
mainly to uplift flow above the guardrails, which pushes the saltation particles to a greater height. 
3.4 Thickness of the sand particle deposition 


The field test shows that the wind-blown sand flux follows a power-law function and the correlation 
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—e— Subgrade with the W-beam central guardrails 
—s— Subgrade without the W-beam central drails 
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Fig. 6 Average saltation height as a function of distance from the W-beam central guardrails. H is the height of the 
subgrade. 


coefficient of the function is between 0.81 and 0.95 (Table 1). The wind-blown sand flux at 1-20, 20— 
40 and 40-80 cm heights accounts for 64.50%—82.31%, 9.87%-17.26% and 7.82%-—18.24% of the 
total wind-blown sand flux in the 0—80 cm sand layer, respectively (Table 1). 


Table 1 Parameters of the wind-blown sand flux structure in the field test 
Percentage of the sand flux to the total flux in the 0-80 cm 


Underlying surface Powerlaw:fit R sand layer (%) Sand transport 
type 0-20 cm 20-40 cm 40-80 cm (ke/(hm)) 
Mobile sand bed y=169.36x 1° 0.81 69.07 15.54 15.38 73.1 
Windward shoulder y=615.6x!35 0.95 82.31 9.87 7.82 17.1 
Leeward shoulder y=253.73x®°®° 0.92 64.50 17.26 18.24 15.3 


To estimate the sand transport in the 0—80 cm layer, we integrate the power-law function of the 
wind-blown sand flux structure. The sand transport is 73.1 kg/(h-m) over the sand bed and 17.1 
kg/(h-m) at the windward shoulder. The sand flux on the leeward shoulder is 15.3 kg/(h-m). The sand 
transported over the sand bed exceeds that transported over the subgrade surface, because the sand 
bed provides a greater source of the sand particles. Some of the sand particles are deposited under the 
influence of the subgrade slope and there is no supply of new sand, so the concentration of the sand is 
less than that on the sand bed. 

Taking the sandstorm occurred on 11 May 2016 at K50 of the Bachu-Shache desert expressway as 
an example (Fig. 7), the inlet and outlet sand mass was 17.1 and 15.3 kg/ (h-m), respectively, and the 
bulk density of the sand particles was about 1.65 g/cm?. The sandstorm lasted about 10 h, and the 
modeled depth of the sand was about 1.1 cm, which was calculated from Equation 3. 


Fig. 7 Photograph of the sand deposition at K50 of the Bachu-Shache desert expressway 
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The field investigation shows that this sandstorm deposited sand in a 0.2—0.3-cm-thick strip and 
that the volume of the deposition depends on both the amount of available sand and the wind power. 
Sand deposition occurs mainly on the leeward side of the guardrails in the overtaking lane. In the 
field test, the distance of the overtaking lane from the W-beam central guardrails was between 65 
and 440 cm. Deposition of the sand on the expressway is consistent with the minimum wind 
velocity and airflow in the wind tunnel test. 


4 Discussion 


For a desert expressway with the W-beam central guardrails, the airflow is affected by the 
combined action of the subgrade and W-beam central guardrails. On the subgrade surface around 
the W-beam central guardrails, the approaching flow is compressed and separates into two flows: 
the bleed flow and the lift flow. The bleed flow passes through the gaps at the base of the W-beam 
central guardrails and is therefore compressed, and its velocity increases slightly. The lift flow 
passes over the top of the W-beam central guardrails, and its velocity also increases. The lift flow is 
separated at the downwind edge of the W-beam central guardrails, leading to a reversed flow in the 
lee of the W-beam central guardrails. The reversed negative flow meets the positive bleed flow 
from the gaps at the base of the W-beam central guardrails, at which point the reversed airflow 
vanishes and the positive wind velocity decreases. Airflow fields behind porous fences are 
complicated by both the bleed flow that passes through gaps in the fence and the displaced airflow 
passing over the fence (Dong et al., 2007). Finally, the eddy zone disappears when the fence 
porosity exceeds 40% (Lee and Kim, 1999). 

Variations in the wind velocity can be used to represent the energy of the airflow, with higher wind 
velocity leading to more wind-blown sand. If sufficient sand is available, then the sand discharge is 
proportional to the third power of the wind velocity (Bagnold, 1941). The subgrade and W-beam 
central guardrails exert a drag force on the wind field, which causes a net loss of momentum in the 
incompressible airflow and therefore a shelter effect. The additional drag force of the W-beam 
central guardrails increases the friction velocity and the roughness element. The energy lost by the 
airflow decreases, so the lower momentum of the aeolian sand leads to disequilibrium in the 
wind-blown sand, and the sand particles settle out in the low-velocity zone (Thomas et al., 2014). 

The wind-blown sand flux is not positively correlated with height, and the maximum mass 
density of the sand and the particle saltation height both increase, similar to the results from the 
Gobi region reported by Zhang et al. (2015). The wind-blown sand flux is affected by the inlet wind 
velocity, the grain size, the fetch length and the characteristics of the underlying surface (Dong et 
al., 2002; Liu and Dong, 2004; Thomas et al., 2014). 

The average saltation height and the decay rate of the wind-blown sand both increase with wind 
velocity and fetch length (Dong, 2004). Sand transportation at the height of 0-2 cm exceeds 50% of 
the sand transportation at the height of 0-10 cm, which means that the near-surface wind-blown 
sand flux is more saturated and there is more active exchange on the sand bed and more sand 
deposits near the W-beam central guardrails (Wu, 1987). The value of the average saltation height 
increases to 2.72-5.17 cm in the range from —7.0H to 20.0H, which shows that the wind-blown 
sand flux on the subgrade surface is not saturated because there is an insufficient supply of the 
sand, indicating the occurrence of saltation. The wind-blown sand is not in active exchange with 
the subgrade surface. The loss of airflow energy near the W-beam central guardrails is the main 
reason for the deposition of the sand. 

In the field test, the structure of the wind-blown sand flux shows that the sand flux activity 
concentrates mainly near the surface in the 0-20 cm layer, which is consistent with the results of 
Cheng et al. (2015). Because of the difference in scale, the characteristics of the wind-blown sand 
flux in the wind tunnel test differ from those observed in the field test. In addition, the lower 
chamber of the collector of the sand sampler in the field was 2.0, 5.0, 7.5 and 9.0 cm deep and the 
cross-sectional area of the chamber was 1.7 cmx1.7 cm. The sparsely spaced collectors and the 
large cross-sectional area of the chamber cannot accurately reflect the near-surface wind-blown 
sand flux. 
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This study investigates the characteristics of the sand deposition on the Bachu-Shache desert 
expressway. The field investigations show that the wind-blown sand deposits mainly on the 
leeward side of the W-beam central guardrails in the overtaking lane. The velocity characteristics 
observed in the wind tunnel test show that the minimum wind velocity and the lowest airflow 
energy occur at 3.5H. The deposition of the wind-blown sand on the desert expressway is consistent 
with the minimum wind velocity and airflow observed in the wind tunnel test. 


5 Conclusions 


In this study, the wind tunnel test shows that the subgrade of desert expressway produces two zones 
of low wind velocity, one on the windward side of the subgrade and the other on the leeward side 
of the subgrade foot. For subgrades with the W-beam central guardrails, a zone of increasing wind 
velocity forms above the W-beam central guardrails and a zone of decreasing wind velocity forms 
in the lee of the W-beam central guardrails. The maximum transport of the sand under the influence 
of the subgrade is concentrated between 3 and 5 cm above the subgrade surface, and the sand 
particles are uplifted because of the uplift airflow and stronger bounce from the subgrade surface 
and W-beam central guardrails. The field test shows that the W-beam central guardrails cause the 
deposition of the wind-blown sand leeward of the guardrails in the overtaking lane, which is 
consistent with the point of the minimum wind velocity in the wind tunnel test. 

This study clarifies the interaction between the W-beam central guardrails and the movement of 
the wind-blown sand, but it does not propose a suitable type of the W-beam central guardrails for 
the desert expressways. Future research should focus on improving the structure of the W-beam 
central guardrails not only to reduce the impact of the wind-blown sand but also to satisfy the 
strength requirements of the traffic safety infrastructure of desert expressways. 
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